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Abstract 

The correlation between the fast cage dynamics and the structural 
relaxation is investigated in a model polymer system. It is shown that 
the cage vibration amplitude, as expressed by the Debye- Waller factor 
(DW), and the relaxation time collapse on a single universal curve 
with a simple analytic form when the temperature, the density, the chain 
length and the monomer-monomer interaction potential are changed. For 
the physical states with the same the coincidence of the mean-square 
displacement, the intermediate scattering function and the non-gaussian 
parameter is observed in a wide time window spanning from the ballistic 
regime to the onset of the Rouse dynamics driven by the chain connectiv- 
ity. The role of the non-gaussian effects is discussed. Keywords: glass 

transition; polymers dyamics; Molecular Dynamics simulation 



1 Introduction 

There is a growing interest in the relation between the fast vibrational dynamics 
and the long-time structural relaxation occurring in viscous systems and super- 
cooled liquids close to their glass transition [Il[l[31[l[51[51[3[Hl[ni[ini[ni[Il[Il 
[11 [m [III 113 [m [H [20]. This resuhed in studies of the vibrational dynamics 
of both glasses [i [13l [H [15l [M [HI and fluid systems [HI H [20] . 

Recently, the universal correlation between the amplitude of the caged dy- 
namics, as expressed by the Debye- Waller factor (DW), and the structural re- 
laxation time Tq has been evidenced by simulations and experiments on several 
physical systems including molecular liquids, polymers and metallic alloys in a 
wide range of fragility [20j. It was found that the shape of the related scaling 
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function is also controlled by the non-gaussianity of the kinetic unit displacement 
[20] . The present paper provides further insight on that issues by presenting 
new numerical results. 



2 Theory 

The glass transition has been pictured as the freezing of a liquid in an Aperiodic 
Crystal Structure (ACS) where the viscous flow is due to activated jumps over 
energy barriers Ai? oc kBTa? / {u^) where a is the displacement to overcome the 
barrier and {u^) is the DW factor of the liquid, i.e. the amplitude of the rattling 
motion within the cage of the nearest neighbours atoms [3] . The ACS picture 
leads to the Hall-Wolynes (HW) equation, Tq,?? oc exp(a^/2(M^)), a relation 
which has been derived by others too [TH [51 [H]. The HW equation relies 
on the condition that Tq exceeds the vibrational time scales. When the HW 
equation is compared with the experiments, one notes strong deviations from 
the predicted hnear dependence between logr^ and [H]. To generahze 

the HW equation, one considers a distribution p{a^) of the square displacement 
to overcome the energy barriers |20j . Following the central limit theorem, the 
gaussian form is adopted, p{a^) = M ■ exp[— (a^ — a^)^/2cr^2] where M is the 
normalization constant, and da'^ are the average and the standard deviation 
respectively. Note that, since Ai? oc a?, this choice corresponds to a gaussian 
distribution of energy barriers. The distribution is taken to be independent 
of the state parameters because the average displacement of the kinetic unit 
within Ta is weakly dependent on Tq, [22]. In contrast, the DW factor depends 
on the state parameters [HI [Z]- If one averages the HW expression over p{a?), 
one obtains the following generalized HW expression (GHW): 

r.=roexp(^^ + ^j. (1) 

Eq. [1] neglects the very weak DW-dependence of tq. If the linear temperature 
dependence of the DW factor is assumed, GHW reduces to other results reported 
for supercooled liquids [23] , polymers |24j or models of glassy relaxation [551 12S] ■ 



3 Model 

A coarse-grained model of a linear freely-jointed polymer is used. Non-bonded 
monomers interact via a generalized Lennard- Jones pair potential Up q{r) with 
Up,,ir) = e{q-p)-^[p{a*/r)i - q{a*/r)P] + U,ut with a* = 21/6^.' The pa- 
rameters p and q control the stiffness of the attractive and the repulsive part, 
respectively. All quantities are in reduced units: length in units of a, temper- 
ature in units of e/ks, and time in units of a^/ m/e, where m and ks are the 
monomer mass and the Boltzmann constant, respectively. The energy unit is 
given by the depth of the potential well e. We also set m = fc^ = 1. The poten- 
tial is cut and shifted to zero by Ucut at r = 2.5. The potential Up.q{r) reduces 
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to the usual Lennard- Jones (LJ) potential by setting p = 6, g = 12 . Bonded 
monomers interact with a potential which is the sum of the FENE (Finitely 
Extendible Nonlinear Elastic) potential and the LJ potential (see ref. [20] for 
further details). This results in a bond length b = 0.97. Samples with N ~ 2000 
monomers were used. Equilibration runs were performed in isothermal-isobaric 
(NPT) or canonical (NTV) ensembles. Data were collected under microcanoni- 
cal conditions by integrating the equations of motion with a reversible multiple 
time steps algorithm, i.e. the r-RESPA algorithm [57]. Physical states with 
different values of the temperature T, the density p, the chain length M and 
the monomer-monomer interaction potential Up^q{r) were studied. See ref. [20j 
for further details. 

4 Results and discussion 

In this section the translational dynamics and the relaxation of the monomeric 
unit are studied. Changing the state parameters (T, p, M, p, q) results in changes 
of both the DW factor and the relaxation time r^. It was found that, when 
different physical states have the same relaxation time, both their translational 
dynamics, as expressed by the mean square displacement (MSD), and their 
relaxation, as expressed by the self part of the intermediate scattering function 
(ISF), are coincident from the ballistic regime up to the onset of the connectivity 
effects (Rouse regime) at times fairly longer than Tq, [50] . The resulting clusters 
of curves are shown in FigHlfor both MSD ( {r^{t)) = A^"^ Ef ((i"i(^) -ri(0))^) 
) and ISF ( Fs{qrnax,t) = N''^ X]f (exp[iqma2;(ri(t) -ri(0))]), qmax refers to the 
maximum of the static structure factor). '[Insert figure 1 about here]' Fig[T]also 
shows the definition of Tq via the equation Fs{qmax,Ta) — 1/e. The existence 
of clusters of physical states with similar dynamics over the wide time range 
from the vibrational regime to the long-time relaxation suggests that the latter 
are correlated. To investigate this issue, it was noted that the location of the 
inflection point of the MSD does not depend on the state point, i.e. it always 
occurs at the same time t — t* \n the log- log plot of MSD [50] . This gives the 
opportunity of a clear-cut definition of the DW factor as {u^) = (r^(t*)). The 
plot of Tq vs (w^) leads to a master curve for the numerical results well described 
by the GHW EqUwhich fits the experimental results over about eighteen orders 
of magnitude [20] . 

We investigated if, in addition to MSD and ISF, other quantities exhibit 
identical time-dependence when evaluated for the clusters of states with identical 
Tq. values shown in FiglTj The results for the non-gaussian parameter (NGF) 
oi2{t) = (3(r'*(t)}/5(r^(t))^) — 1, which quantifies the dynamical heterogeneity 
of the system at a given time t [2S], are shown in Fig[5| '[Insert figure 2 about 
here]' NGF vanishes if the monomer displacement is a gaussian process. Fig|5] 
shows that NGF increases between the end of the ballistic regime (t ~ 0.1) 
and roughly Tq. It is apparent that, within the statistical uncertainty, physical 
states with coincident MSD and ISF (see Fig[T]) have equal NGF up to about 
Tq, too. For longer times the NGF of physical states with equal Tq but different 
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chain length differ also due to the onset of connectivity effects (Rouse regime) 
[291 [30]. 

It has been already shown that the magnitude of the non-gaussian parameter 
is related to the curvature of the GHW Eq[T]when the plot logTc vs {u^)~^ is 
considered [5^. Indeed, the ratio of the quadratic and the linear terms of Eq[T] 
with respect to (u^)^^, TZ = / {^a^ {u^)) increases with the height of NGF 
Oi2max and, if the latter vanishes, TZ does the same [20]. The inset of Fig(2] 
shows the increase of max by increasing Tq . If the non-gaussian effects are 
missing, ISF reduces to F^{qmax,t) = exp{—^q'^^^{r'^)). The first correction to 
F^{q,t) due to the non-gaussian effects depends on NGF and reads [3T| : 

=exp (-Iq^ir^) 

(2) 

'[Insert figure 3 about here]' Fig [3] compares the numerical results for ISF with 
the approximation Fs{qinax,t)- It is seen that, when > 10^, Fs{qrnax,t) 
poorly approximates the numerical results, thus showing that the non-gaussian 
effects on the relaxation are not accounted for by the first correction to F^ {qmax , t) , 
i.e. they are not small. However, in spite of the large deviations of Fs{qmax, t) 
from the exact results, the inset of Figl3] shows that the relative error be- 
tween Ta and the approximated estimate Ta (to be defined by the equation 
Fsiqmax,Ta) = 1/e) is reasonable. 

5 Conclusions 

The correlations between the fast dynamics of the monomers within the cage 
of the first neighbours and the long-time structural relaxation are studied. It 
is shown that physical states with equal Tq exhibit coincident MSD, ISF and 
NGF from the ballistic regime up to the onset of the connectivity effects (Rouse 
regime) at times fairly longer than Tq . The first correction to the gaussian ap- 
proximation of ISF disagrees from the numerical results for Ta ^ 10^. However, 
the relative error between Tq, and the approximated estimate Tq, stays within 
reasonable bounds. 

Acknowledgment 

Discussions with Dr. Luca Larini are acknowledged. Computational resources 
by "Laboratorio per il Calcolo Scientifico" (Physics Department, University 
of Pisa), financial support from MUR within the PRIN project "Aging, fluc- 
tuation and response in out-of-equilibrium glassy systems" and FIRB project 
"Nanopack" are acknowledged. 




4 



References 

1] C.A. Angell, Science 267 1924 (1995). 

2] J. Shao and C.A. Angell, Proc. XVIIth International Congress on Glass, 
Beijing, 1 311 (1995). 

3] R.W. Hall and P.O. Wolynes, J. Chem. Phys. 86 2943 (1987). 

4] J.C. Dyre, N.B. Olsen, and T. Christensen, Phys. Rev. B 53 2171 (1996). 

5] J.C. Dyre and N.B. Olsen, Phys. Rev. E 69 042501 (2004). 

6] L.M. Martinez and C.A. Angell, Nature 410 663 (2001). 

7] K.L. Ngai, .J. Non-Cryst. Solids 275 7 (2000). 

8] K.L. Ngai, Phil. Mag. 84 1341 (2004). 

9] F. Starr ct al., Phys. Rev. Lett. 89 125501 (2002). 



P. Bordat et al., Phys. Rev. Lett. 93 105502 (2004). 

A. Widmer-Cooper and P. Harrowell, Phys. Rev. Lett. 96 185701 (2006). 

U. Buchenau and R. Zorn, Europhys. Lett. 18 523 (1992). 

T. Scopigno et al.. Science 302 849 (2003). 

A.P. Sokolov et al., Phys. Rev. Lett. 7 2062 (1993). 

U. Buchenau and A. Wischnewski, Phys. Rev. B 70 092201 (2004). 

V.N. Novikov and A.P. Sokolov, Nature 431 961 (2004). 

V.N. Novikov, Y. Ding, and A.P. Sokolov, Phys.Rev.E 71 061501 (2005). 

S.N. Yannopoulos and G.P. Johari, Nature 442 E7 (2006). 

J.C. Dyre, Rev. Mod. Phys. 78 953 (2006). 

L. Larini et al.. Nature Physics 4 42 (2008). 

V.N. Novikov and A.P. Sokolov, Phys. Rev. E 67 031507 (2003). 

C.A. AngeU, J.Non-Crystalline Sol. 131-133 13-31 (1991). 

H. Bassler, Phys. Rev. Lett. 58 767 (1987). 

J.D. Ferry, L.D.J. Grandine, and E.R. Fitzgerald, J. Appl. Phys. 24 911 

(1953). 

C. Monthus and J.P. Bouchaud, J. Phys. A: Math. Gen. 29 3847 (1996). 
J.P. Garrhan and D. Chandler, Proc. Natl. Acad. Sci. 100 9710 (2003). 



5 



[27] M.E. Tuckerman and B.J. Berne, J. Chem. Phys. 95 8362 (1991). 

[28] S.C. Glotzer and M. Vogel, Phys. Rev. E 70 061504 (2004). 

[29] G.D. Smith et al., J. Chem. Phys. 114 4285 (2001). 

[30] W. Paul, Chem. Phys. 284 59 (2002). 

[31] B.R.A. Nijboer and A. Rahman, Physica 32 415 (1966). 



6 



Figure 1: (a) MSD time-dependence in selected cases. See ref.piT for details. 
The MSDs are multiplied by indicated factors, (b) corresponding ISF curves. 
Four sets of clustered curves (A through D) show that, if states have equal 
Tq, (marked with star on each curve ), the MSD and ISF curves coincide from 
times fairly longer than Tq, down to the crossover to the ballistic regime at least. 
Adapted from ref. 20] . 



Figure 2: Time-dependence of the the non-gaussian parameter (NGF) for the 
same states of FiglT] The stars denote the time t — Ta- The plot shows that 
for states with equal Tq not only MSD and ISF coincide between t* and Ta 
(see Figd]) but also NGF does the same within the statistical uncertainty. The 
inset shows the increase of the maximum of NGF with Tq . The dashed line is a 
parabolic guide for the eyes. 



Figure 3: Comparison between the numerical ISF (continuous lines) of Fig[T] 
and the first correction to the gaussian approximation Fs {qmax , t) (dashed lines) . 
The inset shows the relative error between Tq, and the estimate Tq, as drawn 

from Fs{qmax,Ta) = 1/e. 
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FIGURE 2 
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FIGURE 3 
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